GLOBAL RESULTS FOR LINEAR WAVES ON EXPANDING 
SCHWARZSCHILD DE SITTER COSMOLOGIES 



VOLKER SCHLUE 

Abstract. In this global study of solutions to the linear wave equation on 
C "3 l Schwarzschild de Sitter spacetimes we attend to the cosmological region of 

O^l ■ spacetime which is bounded in the past by cosmological horizons and to the 

future by a spacelike hypersurface at infinity. We prove an energy estimate 
capturing the expansion of that region which combined with earlier results for 
the static region yields a global boundedness result for linear waves. It asserts 
1 that a general finite energy solution to the global initial value problem has a 

^v^j 1 limit on the future boundary at infinity that can be viewed as a function on 

the standard cylinder with finite energy, and that moreover any decay along 
the cosmological horizon is inherited along the future boundary. In particular, 
O ■ we exhibit an explicit nonvanishing quantity on the future boundary of the 

spacetime consistent with our expectations for the nonlinear stability problem. 
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1. Introduction 

The primary impetus for the study of the linear wave equation on black hole 
spacetimes is that some aspects of the stability or instability of dynamical black 
holes are already foreshadowed in the global behaviour of linear waves [DR11, 
lArellj . In this paper, we affirm the global linear stability of the simplest cosmo- 
logical black holes, namely spherically symmetric subextremal vacuum black holes 
with positive cosmological constant. While some important aspects of this problem 
pertaining to the linear stability of the static domains have already been resolved 
in |BH08| IDR07| IMSV08| |Dyall| IVaslOj , we focus here on the cosmological region 
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Figure 1. Global geometry of Schwarzschild de Sitter. The ex- 
panding region is the domain bounded in the past by the cos- 
mological horizons C + and C+ and to the future by the spacelike 
hypersurface E + at infinity. H + and H + denote the event horizons 
of the adjacent black hole interiors. 

to the future of the static domains undergoing accelerated expansion. We identify 
a global redshift effect as the relevant stability mechanism which combined with 
previous work for the static region yields a global linear stability result. 

The cosmological black hole spacetimes underlying the present study are pre- 
cisely the members of the Schwarzschild de Sitter family in the subextremal range. 
(Slowly rotating Kerr de Sitter cosmologies will be treated as a perturbation in a 
subsequent paper.) We may depict these spacetimes as an infinite chain of black 
hole interiors, static black hole exteriors, and cosmological regions separated by 
event horizons and cosmological horizons respectively JGH77] , The distinctive fea- 
ture in comparison to the case of vanishing cosmological constant is here that future 
null infinity is a spacelike hypersurface E + , along with the existence of cosmological 
horizons C + , C + as future boundaries of the static regions in the outgoing direction. 
The domain bounded in the past by the cosmological horizons C + U C + and to the 
future by E + has the property that the area of each sphere is strictly increasing 
towards the future and shall thus be referred to as the expanding region 1Z + ] see 
figure [TJ 

The expanding region is foliated by spacelike hypersurfaces of constant area 
radius. We exhibit in this paper an energy of linear waves which is monotone with 
respect to this foliation. Moreover, we show that any solution to the linear wave 
equation with initially finite energy on a Cauchy hypersurface E lying in the past of 
E + (such that E + is in the domain of dependence of E, see figure [T]) attains a limit 
on E + which can be viewed as a function on the standard cylinder lx§ 2 with finite 
energy. In addition, it will be shown that as a consequence of the global redshift 
effect, and the local redshift effect on the horizon, polynomial or exponential decay 
along the cosmological horizon C + will be inherited along the future boundary E + . 
The ideas of the proof originate in the author's treatment of linear waves on de 
Sitter cosmologies |Schl2) referred to in the epilogue of this paper. 

While the resolution of the black hole stability conjecture has remained elusive 
for many years |Daf09] . the motivation for the linear theory developed in this pa- 
per is largely drawn from our expectation that the nonlinear stability problem is 
within reach in the context of expanding spacetimes. The favorable role of a pos- 
itive cosmological constant, which allows for expanding cosmological solutions to 
the Einstein equations, was first recognized by Friedrich in his proof of the stabil- 
ity of the de Sitter solution |Fri86j . While the conformal method has since then 
been extended to scalar field models |Rin08[ I AMV] . the stabilizing effect of the 
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cosmological constant may be seen most clearly manifested in the recent works of 
[RS09[ [Spel2| on perfect fluids, (see also |LVllj for the pure radiation case). 

In view of these results for the nonlinear field equations with positive cosmolog- 
ical constant and the linear theory presented in this paper, we expect the charac- 
teristic initial value problem for perturbations of Schwarzschild de Sitter data on 
C + U C + to be more tractable than the major open problem of establishing the 
nonlinear stability of the past of Ti + U C + ; c.f. figured] 

2. BOUNDEDNESS AND DECAY RESULTS 

2.1. Geometry of the Schwarzschild de Sitter family. The black hole space- 
times of interest in this paper are the members of the Schwarzschild de Sitter family 
(-M(A,m)j fl 1 ) f° r which A > and < 3m\/A < 1 (where A denotes the cosmologi- 
cal constant and m the mass of the black hole). An important role in our analysis 
and in the depiction of the causal geometry of this spherically symmetric spacetime 
family is played by the function r on Q = .M/SO(3) which takes the value of the 
area radius on a sphere of symmetry. In other words, r is constant on the orbits 
of the SO(3) subgroup of the isometry group of the spacetime (■M(A, m ) ! S'). In the 
given parameter range the spacetime (M(h,m)i d) is partitioned into black/white 
hole interiors B, static exteriors S and cosmological regions TZ in the sense that 



where a sphere is in B, S or TZ according as to whether its area radius is in the 
interval < r < r% , r% < r < rc or r > rc respectively; here 2m < r% < 3m < rc 
are the positive roots of a cubic polynomial with coefficients in A and m (for explicit 
expressions see [LR77] ) which mark the event, Ti,T-L, and cosmological horizons, C, 
C, respectively. 

In the following we fix A > 0, < 3m < 1/vA and pick a domain TZ; as discussed 
in [LR77, Schl2] — and included for completeness in Section llOl — the cosmological 
region TZ and its adjacent static domains S can be covered with a single double null 
coordinate system (u, v) (i.e. u, v are functions on Q increasing towards the future 
such that the level sets of u are outgoing, and the levels sets of v are ingoing null 
lines) which identifies the cosmological horizons C, C, with the null hypersurfaccs 
v = 0, u = 0, respectively (see figure [2]). We refer to the component of TZ to the 
future of the static regions S as the expanding region TZ + which in the above chart 
corresponds to < uv < 1. The hyperbola uv = 1 represents the future timelike 
boundary of TZ + : a spacelike cylinder composed of spheres of infinite radius. We 
shall also denote by C+ = C n {u > 0}, C+ = Cl~l {v > 0}, the components of C, C to 
the future of 5, respectively, so that the past boundary of TZ + is precisely C + U C + . 

In the expanding region the area radius r > rc is strictly increasing along any 
future-directed timelike curve. We shall thus view r as a time function for the 
expanding region and obtain a foliation of TZ + by the level sets E r of r : TZ + — > 
(rc, oo ) with the lapse function 



where = —g^ v d v r denotes the gradient vectorfield of r. The metric g decom- 
posed in this way reads 



(2.1) 



M = ...UBUSUTZUSUBU... 



(2.2) 



1 



1 




(2.3) 



g = -4> 2 dr 2 + g, t 
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Figure 2. Causal geometry of the cosmological region 1Z = {0 < 
uv < 1} and the adjacent static exteriors S = {uv < 0}. The 
expanding region 1Z + is the component of 1Z to the future of the 
cosmological horizons C — {v — 0} and C = {u = 0}, and the 
future boundary E + of 1Z + is represented by the hyperbola uv = 1. 



where g r is the induced metric on E r ; in fact we can choose a function t : 1Z + — » 
(—00,00) which is constant on the lines u — av, a € (0, 00), such that 

/Ar 2 2m \ , , 9 o 

(2.4) g r = (— + — -ljdt 2 + r 2 7, 

o 

where 7 is the standard metric on the unit sphere. The level sets S r are thus 

o o 

conformal to the standard cylinder (R x § 2 ,3= dt 2 + 7), and we note in particular 
for the volume form 

(2.5) <p dfig r = r 2 d/io . 

The normal to E r is given by n = 4>V , while we denote by 

(2.6) ^- = <fv = cj ) n 

the vectorfield colinear to V with the property that d r r = 1; (in the above (t,r) 
coordinates d r is indeed the coordinate vectorfield). 

Furthermore there exists a global Killing vectorfield T which coincides with the 
stationary vectorfield in the static regions (where it is timelike), and with the co- 
ordinate vectorfield in the expanding region (where it is spacelike); its integral 
curves in the static region correspond to physical time translations, and its prop- 
erties on the cosmological horizons (where it is null) are intimately linked to the 
positivity of the surface gravity (see Section [32]) . Note that T and the generators 
of the spherical isometry : i — 1, 2, 3 (which are discussed in Appendix fAI are 
orthogonal to n; the tangent space to E r is thus spanned by Killing vectorfields at 
each point. 

2.2. Linear wave equation. This paper is devoted to the global study of solutions 
to the linear wave equation 

(2.7) D g tp = 

on the family of black hole spacetimes (.M(A,m), g) described in Section |2~T1 Our 
results apply to a class of finite energy solutions for the Cauchy problem associated 
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to (|2.7j) . Let us recall for this purpose briefly the standard energy current J x 
associated to the multiplier vectorfield X, a 1-form defined by 

(2.8) J x [ip] -N = T(X,N), 

where T denotes the standard energy momentum tensor associated to (|2.7p (for 
an introductory discussion of Lagrangian field theories in a curved spacetime see 
Chr08 ). Given a spacelike hypersurface (£,<?) with normal n we can think of 

(2.9) / J n [4i] •nd/% > 

JT, 

as the total energy flux of the wave tp through S. Our assumptions on the solutions 
of (|2.7[) considered in this paper are stated in terms of a finiteness condition on the 
level of initial data for an expression of the form (|2 .9[) . 

2.3. Integral estimates. Our results provide integral bounds for the energy den- 
sity of linear waves on spacelike noncompact hypersurfaces in the expanding region, 
in particular on the future boundary, in terms of the energy on general spacelike 
hypersurfaces, and thus emphasize the global nature of our approach. 

In Section H3] we demonstrate the presence of a "global redshift" which leads us 
to our key result for the expanding region. 

Proposition 2.1. Let ip be a solution to (|2.7p then for all r 2 > r\ > tq we have 




Notation. Here denotes the projection of the covariant derivative V of g to the 
sphere of area radius r at this point, and we simply write |-| for the induced norm 
on this sphere. 

It is essentially due to the local redshift effect of the cosmological horizons that 
Prop. [2TT1 can be turned into a global result. We should regard the positivity of the 
surface gravity of the cosmological horizons (see Section I3.2[) as the origin of this 
effect which is here manifestly exploited in two ways. 

Firstly in a construction of a redshift vectorfield N (we refer to the epilogue 
in |DR08] for a general discussion) which gives rise to an energy identity on the 
spacetime domain bounded in the future by S ri and in the past by C + U C + ; see 
Section E3~2"l By choosing r\ > rc sufficiently small, we obtain an estimate for the 
right hand side in Prop. [2~T1 in terms of the corresponding energy on the cosmological 
horizons. In fact, we have: 

Proposition 2.2. Let £ C TZ + be a spacelike hypersurface with normal n crossing 
the cosmological horizons to the future of C + l~l C + , and let us denote by C$ = 
C + D J + (E) the segment ofC + to the future o/S. There exists a timelike vectorfield 
N (which is invariant under the 1-parameter group of isometries <pt generated by 
T) in a neighborhood of of the form M = {r^ < r < r^} fl J + (S) (where 
Tq < rc < r^ are constants that only depend on S, A and m ) with the property 
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Figure 3. Notation for hypersurfaces in Propositions 12. 14 l 2 . 31 



that for all solutions ip to the wave equation (|2.7[) , 



E ri nJ+(S) 



^> 2 + ^<f> 2 + lwl 2 W^ 



<C(n) / *J N [iP] + C( ri ) [ J N [^}-ndfig. 

JC+ JT,' 



whenever rc < r\ < (provided r\ is chosen small enough so that E ri PI £ 7^ ) 
where C is a constant that only depends on r\ and £' = J _ (S ri ) (1 £ denotes the 
spacelike hypersurface £ truncated by S ri . 

Notation. Given a 1-form J we denote its dual with respect to the volume form by 
* J. The notation g refers to the first fundamental form of S. 

Secondly in the work of [DR07] which is occupied with the construction of suit- 
able currents for the static region. We cite here a boundedness result which yields 
control on the nondegenerate energy on the cosmological horizon in terms of a nat- 
ural energy on a spacelike Cauchy hypersurface and thus closes our estimate as 
desired. 

Proposition 2.3 (Prop. 10.3.2, DR07 ). Let £ c S be a spacelike hypersurface 
with normal n in the static region crossing the horizons to the future of the bifur- 
cation spheres C (1C, H fl H, and denote by Cq = C + fl J + (S). We have for all 
solutions ip to (|2.7p . 

*J N m<C f J^-nd^, 
c+ Jt. 

where C is a constant that only depends on E ; A and m, whenever the right hand 
side is finite. 



Remark 2.4. Here N is precisely the vectorfield referred to in Prop. 12.21 While 
DR07^ relies on an extension of N into the static region, we require for our purposes 
an extension into the expanding region; however, on the cosmological horizon C + 
the two constructions coincide. 



Remark 2.5. We point out that a remarkably short proof of Prop. 12.31 is given in 
DR08 (see in particular Corollary 7.1 therein) separating the boundedness from 
the decay argument of [DR07]- The implications of related decay results for the 
static region |DR07[ IEH081 IMSV081 |DyalO[ |Dyall| are discussed in Remark 

Let us consider global solutions to the Cauchy problem (|2.7|) with initial data 
prescribed on a spacelike hypersurface S C S U 1Z + as depicted in figure [1] more 
precisely we say: 



LINEAR WAVES ON SCHWARZSCHILD DE SITTER 



7 



Definition 2.6. A spacelike hypersurface E C M.(\,<m) is called a E + -Cauchy 
hypersurface if the future boundary E + of a chosen expanding region 1Z + C A^(A,m) 
is in the domain of dependence of E C J~(E + ), and E C TZ + US crosses the 
horizons H, C, C, and H to the future of the bifurcation spheres HDH, C DC, and 
HnH, respectively. 

As a consequence of the global and local redshift effect on and near the cosmolog- 
ical horizons — namely Propositions 12.1112.21 and 12.31 — the energy on E r n </ + (E) 
is hnite for any r > rc provided it is initally finite on E. We may thus take the 
limit r — > oo in the global redshift Proposition 12.11 which yields the finiteness of an 
explicit limiting quantity on E + . In fact, we infer from (|2.5|) that 

is the rescaled asymptotic energy of a solution to the Cauchy problem viewed as a 

o 

function on the standard cylinder; (here ft denotes the covariant derivative on the 
standard sphere (§ 2 ,7)). 

Theorem 2.7. Let E + be the future boundary of a chosen expanding region 1Z + C 
■M(A.m) of Schwarzschild de Sitter (with < 3m-\/A < 1, A > 0) endowed with the 

o 

standard metric g of the cylinder lx§ 2 , and let E be a E + -Cauchy hypersurface 
with normal n. Then all solutions ip to the wave equation (|2.7[) with initial finite 
energy on E, 

(2.11) D[il>]= [ J n [iP]-nd(ig <oo, VlseH^E), 

are globally bounded on 1Z + in the energy defined by Prop. \2.1[ and have a limit on 

o 

E + in H (K x § 2 ). Moreover, the limit as a function on 1 x § 2 satisfies 

(2.12) / |v 7/»| 2 d/i= < C(A,m,'E)D[ij}} , 
Js+ 9 

o 

where C is a constant that only depends on A, m, and E, (and V denotes the 
standard gradient on the cylinder K x S> 2 ). 

Remark 2.8. Our result Thm. 12.71 can be read as an implicit version of a decay 
statement for the "natural" energy in this problem. For we establish the bounded- 
ness of a limiting rescaled quantity which corresponds to the decay of the induced 
geometric quantities arising in this setting; note e.g. on the level of energy densities 
that 

o 

\J7ip\ = r\pil>\ , 

i.e. the angular derivatives on the rescaled standard sphere differ from the induced 
angular derivatives on a sphere on E r by a factor of r. 

Moreover the result is in agreement with our expectation for the nonlinear sta- 
bility problem. Although we expect to recover the same global causal geometry, 
the dynamical development of a perturbation of Schwarzschild de Sitter initial data 
is not expected to settle down to the exact geometry of the expanding region a 
Schwarzschild de Sitter solution. This is captured by a nonvanishing bounded 
quantity on the future boundary. 
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Note that we have dropped the normal derivative in the limit (|2.10p . If we viewed 
the asymptotic quantity recovered in Thm. 12.71 as initial data for the corresponding 
backward problem, then the rescaled asymptotic normal derivative of Prop. 12.11 
would be part of the data at infinity. 

This result should be viewed as the most general global energy bound in the 
class of finite energy solutions to (|2.7p . While in Thm. 12.71 we only require the 
initial energy to be finite, it is known that under suitable assumptions on the 
higher order energies (i.e. under stronger regularity assumptions on the initial data) 
the solutions to (|2.7|) will in fact decay along the cosmological horizons. It is an 
immediate consequence of our approach (see also Section I4.4[) that our estimates 
can be localised to show that any decay along the horizons is in fact inherited along 
the future boundary. 

Notation. Let r be a function on C + which is constant on the spheres of symmetry 
such that t = T > on a chosen sphere S C C + to the future of C + n C + and 



We denote by C+ the segment of C + lying to the future of the sphere with value 
r = n, and by C Tl the outgoing null hypcrsurface from the sphere r — t\\ see 
figure 21 Moreover let us denote bjQ 

s+ = j+(c+uc T )ns + 

the future boundary of the spacetime domain to the future of U C T . 



Corollary 2.9. Assume ip is a solution to (|2.7[) which decays polynomially along 
C + in the sense that for some fixed integer k 6 N there exists a constant C'k > 
such that 

(2.13) / *J W M<% (r>T ). 

Jet T 

Then there exists a constant Ck < C < oo such that also 

(2.14) f |v^| 2 d M o<-^ (r>r ). 



We use standard notation for causal sets in Lorentzian geometry, see e.g. IHE73| . 
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Moreover, if the decay along C + is assumed to be exponential, i.e. if for some ft > 
there exists a constant Cp > such that 

(2.15) / *J N m<C e-^ (r>r ), 



then also the decay along the future boundary S + is exponential, with a rate however 
that is not faster than dictated by the local redshift effect of the cosmological horizon; 
in fact then there exists a constant Cp < C < oo such that 



(2.16) / |v^| 2 d/io <Ce-« T (r>r ), 

where b is a constant that only depends on A, and m; (here for simplicity ft ^= b). 



Remark 2.10. The assumptions of Cor. 12.91 are known to be satisfied under higher 
regularity conditions on the initial data |DR071 [I)yall| . (The constants C^, Cp, 
then depend on higher order norms of the initial data.) 

More precisely, it is shown in |DR07j that (|2.13[) is satisfied with k = m + 1 (for 
any m € N) for any solution to (|2.7[> in the domain of dependence of a E + -Cauchy 
hypersurface £ with normal n (where r = To on £ (1 C + ) provided 

(2.17) f V • n dfig < oo , 

is i 

\a\<m 

where a is a multi-index of order \a\ = ot\ + ot2 + 03 and S1 Q = fl" 1 ^ 2 ^ 3 . 

Moreover, the geometric scattering theory initiated in }SZ97j confirmed that the 
exponential decay assumption (|2.13[) can be satisfied. While [BH08, MSV08 showed 
exponential decay in the class of smooth initial data, Dyatlov |DyalO[ |Dyall| es- 
tablished an exponential decay estimate for solutions to (|2.7[) in the static domain 
holding up to and including the cosmological horizon, thus replacing the polynomial 
decay rate on the right hand side of (|2.13|) by an exponential bound (|2.15[) under 
suitable regularity assumptions on the initial data. Note that the decay results of 
|Dyall hold in the wider class of slowly rotating Kerr de Sitter black holes; see 



also [VaslOj 

2.4. Pointwise estimates. The foliation (S r ) is not only a natural choice for the 
energy estimates in this paper but also particularly convenient from the point of 
view of pointwise estimates. Note that the tangent space to E r (r > re), and 
S + is spanned by Killing vectorfields, namely T and Q.u\ : i = 1,...,3. Since 
the commutations with these vectorfields are controlled in L 2 by Thm. 12.71 we can 
immediately pass to pointwise estimates using the classical Sobolev embedding. 
Indeed, we can alternatively write the result (|2. 12|) using the coercivity equality on 
the sphere (see Appendix [X} as 

r 3 
(2.18) / {(r.^) 2 + ^(%^) 2 }d^<c^]. 

Corollary 2.11. Let £ be a -Cauchy hypersurface as in Thm. \2.7\ and let ip 
be a solution to (|2.7| which in addition to (|2. 1 1|) satisfies on E the higher order 
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energy condition 



3 3 



(2.i9) dm = D[i>] + 53 D[n (i) ^] + 53 D[n w n w) ^] 




+ Dpv] + 53 D[n w T^] + 53 D[n w n W) r^] < oo . 



z— 1 — 1 



TTien we ftave </ie pointwise bound 



(2.20) supl^l < supl^l +C(A,to,£)a/A#] 



on i/ie future boundary T, + , and the pointwise estimates 



(2.21) sup (^) 2 + |X7 V| 2 < C(A, m, E) A#] , 



where C is a constant that only depends on £, A, and m. 

Remark 2.12. If the higher order energies associated to the solutions fl^tp, Qr^Qfj^tp, 
and Ti/>, fi^T^i, f2( i )fi(. ; -)T'!/> are assumed to decay along the cosmological horizon, 
similarly to Cor. I2.9( with either a polynomial rate according to (|2.13j) or exponen- 
tially according to (|2.15|) . then so does the quantity 



along E + , according to (|2.14l) or (|2.16[) respectively. Moreover, pointwise decay for 
tp to a constant follows from (|2.22[) by integration on £+. 

Our result states in particular that solutions to the linear wave equation on 
Schwarzschild de Sitter spacetimes have a limit on the future boundary that can 
be viewed as a function on the standard cylinder lx§ 2 , and moreover, that the 
rescaled induced derivatives have a limit or equivalently that the geometrically in- 
duced derivatives of these solutions decay to the bounded derivatives of the limiting 
function. We have seen that as a consequence of the global redshift effect quanti- 
tative decay that is present along the horizon translates into a quantitative decay 
along the future boundary. 

Remark 2.13. In this paper we have restricted ourselves to spherically symmetric 
spacetimes which has the advantage that the angular momentum operators 
generate exact isometries. In a subsequent addendum to this paper we will treat 
slowly rotating Kerr de Sitter cosmologies as a perturbation of the Schwarzschild 
de Sitter scenario, in fact of the de Sitter cosmology as discussed in Section [5j In 
view of the manifest stability of the redshift construction given in Section 14.11 and 
a smallness condition on the specific angular momentum, similar global results can 
be deduced in this wider class of cosmological black hole spacetimes. 

3. Global geometry of the Schwarzschild de Sitter spacetime 

The Schwarzschild de Sitter spacetimes (M (A,m) > ff) are a family of 3+1-dimensional 
Lorentzian manifolds, which in the range A > 0, < 3m < 1/S/A constitute our 
simplest model for an asymptotically flat black hole spacetime in an expanding 



(2.22) 
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universe; (m — is de Sitter spacetime, our simplest cosmological model of an ex- 
panding universe |NB09j ). They are unique in the sense of Birkhoff [Hcu96 ( DR08 
as spherically symmetric solutions to the vacuum Einstein equations 

(3.1) - -g^R + Ag„„ = 

with positive cosmological constant A > 0. While their global causal geometry is 
well known since |GH77( ILR77] , a few remarks are in order concerning the coordi- 
nates used in this paper. 

3.1. Cosmological spacetimes under spherical symmetry. In spherical sym- 
metry the topology of the spacetime manifold is necessarily that of Q x SO (3) where 
Q is a 1 + 1-dimensional Lorentzian manifold. In a double null foliation of Q by 
the level sets of functions u, v : Q — > R the metric g assumes the form 

(3.2) g = -Q 2 dudv + r 2 7, 

o 

where 7 is the standard metric of the unit sphere § 2 , and we can think of (|3.1j) as 
partial differential equations for the area radius r and the conformal factor Q as 
functions of u, and v. Indeed, it is easily deduced from (I3.1[) that the area radius r 
satisfies the Hessian equations 

(3 3a) ^_l^^: = o 

du 2 SI du du 

d 2 r ldrdr_ il 2 n 2 ^ 

du dv r du dv 4r 4 

d 2 r 2 dfl dr _ 

(3 ' 3c) 

We observe that the mass function m : Q — > K implicitly defined by (and motivated 
in |Chr95Q 

2m Ar 2 4 dr dr 
( ' 1 ~~ 3~ ^~Wdu~dv' 1 

is constant by virtue of (|3.3[) . and precisely the quantity that parametrizes the 
Schwarzschild de Sitter family for any fixed A > 0. It is useful to introduce the 
"tortoise coordinate" 

(3.5) ' ' 



i 2m_ Ar 2 

r 3 



which satisfies again as a consequence of (|3.3[) the simple partial differential equa- 
tion 

3 6 d<2r * 

du dv 

In general, the dependence of the area radius r on the null coordinates u, v, is thus 
of the form 

(3-7) r* =h{u) + h{v), 

where we are free to choose the real valued functions fx, and /2- We may think 
of the different charts on (A4(A, m ), g) to be obtained with suitable choices of the 
functions /i, /2, and of the constants for the indefinite integral (|3.5j) . While in the 
case m = the manifold can be covered with a single chart, this is not possible 



12 



VOLKER SCHLUE 



for m > 0. Here we discuss a chart that covers the cosmological regions and the 
adjacent static domains. 

3.2. Double null coordinates for the cosmological and static regions. The 

horizons of Schwarzschild de Sitter are null hypersurfaces consisting of spheres of 
constant area radius, and thus null lines on Q on which the right hand side of Q3.4p 
vanishes. The polynomial on the left hand side of Q3.4|) has three distinct real roots 
rc, r%, and rc provided 

(3.8) < 3m < -= , 

VA 

satisfying (explicit expressions may be found in [LR77] ) 

(3.9) rc < < 2m< r-u < 3m < r c ■ 
We note 

(3.10) r -2m-^- = -~(r - r n )(r - r c )(r + \fc\) , 
with 

(3.11a) -\rc\ +r c + r n = 

3 

(3.11b) r c \rt\ + r n \rc\ ~ r-urc = 

(3-llc) r H r c \rc\ = — , 

and by decomposition into partial fractions 

i 3 r n i 



(3.12) 



k{r c -r u ){r H + \rc-\)r - r n 
3 r c 1 



A (r c - r H )(r c + |r c |) r - r c 

, 3 Ifcl 



^ (\rc~\ + r<H){\rc\ + r C ) r + \rc-\ ' 

The chart covering the cosmological horizons r = rc is now obtained by choosing 
HI to be centred at r = 3m, 



'• <<>= / , J A^ , 



3 m 



3 



and setting 



(3.14a) A(x) = h(x) = " C ^ log^- 

A (r c - rn)(rc + \rc\) A 

where A is the constant 

(3.15) A 2 = (r c - 3m)(3m - r n y~% (3m + |fc|) _i ^ . 

Indeed, the integration of (13. 13)) using p,12j) yields in view of p,14j) the following 
relation between the chosen null coordinates (u, v) and the radius function r: 

r - r c 



(3.16) uv 



(r — r^i) r c r H + \ r c\ (r + Ire ) rc r d+ r « 
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We recover precisely the picture of figure [2] where the cosmological horizons r = rc 
are located at u = and v = 0, while the surfaces of constant r ^ rc are spacelike 
hyperbolas in the wii-plane for r > rc and timelike hyperbolas for r < rc- Moreover 
the future timelike boundary r = oo is identified with the spacelike hyperbola 
uv = 1. 
Since 

(3 17a) dr * - 1 dr - 3 rc 1 

<9u 1 - 22i _ Atl g u A (r c - r n )(rc + |rc|) « 



(3.17b) 



dr* 1 dr 3 rc 1 



dv i-2m.-Ag.dv ^ {r c ~ r H ){r c + \rc\) v 



we can solve (|3.4[1 for fl 2 to obtain: 
4 3 r 2 



(3.18) Q 



2 



r A (rc - r-u) 2 (r c + \r c \) 2 

x (r - r-Hj rc 7 w + F ^ (r + |fc|) T " c ^ r c\+ r n 

The (nondegenerate) metric g on the chart that covers the region r-u < r < oo and 
extends across the cosmological horizon r = rc thus takes in double null coordinates 
(u, v) the form 

4 3 r 2 

(3.19) q = — t -^t : f-rrX 

r A (rc - r H ) (rc + \rc\) 

*{r-r u ) r c -H + i'-ei (r + |fc|) " c |rcl+ ''" dudw + r 2 7 
where r is a function of (u, u) implicity given by (|3. 161) . 

3.3. Constant area radius foliation of the expanding region. We have seen 
that the expanding region 



(3.20) 1l + = {(u,v) 



< uv < 1 



can be foliated in a geometrically natural way by spacelike hypersurfaces S r of 
constant area radius r, 

oo 

(3.21) 7e+ = (J S r . 

<r<oo 

Here the area radius serves as a time function on 7?. + C M.(\ m \ (r is strictly increas- 
ing along any future-directed timelike curve in 1Z + ). We obtain the decomposition 
of the metric g with respect to the foliation (|3.2ip with the help of the gradient 
vectorficld 

(3.22) V< t = -gTd v r. 
Indeed, since 

1 A (rc - r H )(r c + \rc\) d ^ d 
rc 

we find for the lapse function 



we fine 
(3.24) 
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and the metric takes the form ()2.3[) . namely 
(3.25) g = -4> 2 dr 2 +g r . 

Moreover, we have (|2.6I) . The explicit form of the induced metric g r on S r is readily 
inferred from (|3.19l) upon the introduction of the coordinate 

(3-26) t = -Tl i=S T log| — | . 

A (r c + \r c \)(r c - r H ) 61 u l 

We can then read off (|2 .4[) from the explicit expression for the metric g in (t, r) 
coordinates, which reads 

1 , 9 /Ar 2 2m \ , , 9 o 

( 3 - 27 ) 5 = - Ag.an dr 2 + (— + V " +r 7 ■ 

3 r 

Note here that it follows from (|3.16[) and p. 261) that 

d _ 1 A (r c + |r^K_r c - r H ) ^ 

rc 

and (which we have already made use of in (|3.23p ) 



/ooon 9 1 ^ (rc + \r C \)(rc - r n ) ( d d\ 

(3 - 28) ai = 23 7c \ u d~u ~ v d-v) 



(3.29a) Or = lA (r e + lr e l)(r e -rH) tf 

K ' du 4 3 r c 

ttOQU <9r lA (r c + |rc|)(rc-r^) a 

3.29b — = SI u. 

ov 4 3 rc 

3.4. Surface gravity of the cosmological horizon. In due course the positivity 
of the surface gravity of the cosmological horizons will be of crucial importance. 
While this fact in itself is well-known |GH77) . we take its statement as an oppor- 
tunity to introduce a vectorfield Y which is the starting point for all local redshift 
multiplier constructions (see |DR08j ) . 

The observation is that the vectorfield (13.28[) extends to a global vectorfield that 
characterizes the cosmological horizon as a Killing horizon with positive surface 
gravity. 

Lemma 3.1 (Positive surface gravity of the cosmological horizons). The vectorfield 



is globally Killing, i.e. 



(3.31) ^=^C T g = 0, 
and satisfies 

(3.32) V T T = k c T on C + , 
where 

„ 1 A {r c - r H )(r c + \rc\) v 

(3.33) k c = - — > 

2 3 r c 

is the surface gravity of the cosmological horizons. 
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(3.34) 



The result is obtained with the help of the vectorfield 

2 3 



Y 



|r dv 

ov 



which is conjugate to T along the horizon (i.e. Y is null, orthogonal to the tangent 
space of the spheres of symmetry, and verifies the following normalization): 



(3.35) 

Indeed, by (|3~3l"]l . 
(3.36) g(V T T,Y) 



9(T,Y) 



= -2. 



c+ 



= -g{V Y T,T) 



c+ 



--Y-g(T,T) 



c+ 



d / 2m Ar 2 



dr V r 3 
Alternatively kc is characterized by 

(3.37) VyT = V T Y = -k c Y on C+ ; 

note that this in particular implies that Y is Lie transported by T along the horizon: 



(3.38) 



[T,Y] 



= 0. 



c+ 



Remark 3.2. The vectorfield Y" as defined above forms the basis of an extension 
away from the horizon which captures the redshift effect of horizons with positive 
surface gravity. It was first introduced in the context of the Schwarzschild black 
hole spacetime in |DR09i IDR08] . 

4. Local and global redshift: 
Proofs of the boundedness and decay statements 

We develop in this section the global energy estimates for solutions to the wave 
equation on Schwarzschild dc Sitter spacetimcs. 

The crucial estimate for our analysis is proven in Section |4. II where we present 
a vectorfield that captures that global redshift effect of the expanding region. In 
Section 14.21 the local redshift effect of the cosmological horizons is exploited to 
control the energy in its immediate vicinity. We obtain a global estimate with the 
help of a boundedness result for the static region from |DR07] as discussed above. 

4.1. Global redshift effect. It is the purpose of this section to show that the 
vectorfield 



(4.1) 



», 1 9 

M = -5T 
r or 



captures the global redshift effect in the expanding region of Schwarzschild de Sitter. 



The divergence theorem applied to the energy current 
multiplier M, J M , on 

(4.2) n;i= |J E r , 

ri <r<T2 

combined with the coarea formula applied to the foliation 



associated to the 



(4.3) 



J M ■ ndfj, 



g r2 



dr 



V-J M d/Jg = 



reads 



J M ■nd/fc 
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for any ri > n > rc- 

The redshift property of M is now manifested in the fact that for any solution 
ip to the wave equation (|2.7I) it holds 

(4.4) cj> V • J M [ip] > - J M [ip]-n, 

r 

which is the subject of the proposition below. We are thus in the situation that 
tM i,, , f 2 A J- f tM ^ f jM 



(4.5) / J M -ndfig r + dr- J M -ndfj,g r < / J M -nd^ r 
which implies by a Gronwall inequality: 

(4.6) r 2 / J M ■ ndfe < r x f J M ■ ndte (r a > r x > r c ) . 



9r 2 - ' 1 / " "WA^ 



This is the content of Proposition 12.11 

We begin with proving (|4.4|) . Recall here the conservation laws for solutions to 
the wave equation, 

(4.7) V-T = 0, 
namely 

(4.8) t m „ [if>] = d v <4> - l - 9tlv a a ^d a i> 



is conserved provided ijj is a solution to (|2.7I) . Therefore we obtain for the energy 
current associated to M, 

(4-9) ^[^T^p", 
that its divergence is given by 

(4.10) V • J M = V^jf = ^VT^fofr] = K M [^\ , 

where = ^£m5 is the deformation tensor of M. 

Proposition 4.1. Lei J M be the current associated to the multiplier (|4.1[) . and n 

the normal to S r (r > re)- TTierc /or any solution tp to the wave equation (|2.7|l me 

0V- J M M > i J M [0]-n 
r 

on TZ+. 

Proof. It is equivalent to establish 

r z Z 

Recall from (1X27)) that 

3 r 

<7rr 



A (r - r«)(r - r c )(r + |rc|) 

ff« = ~(r - r«)(r - r c )(r + |fc"|) 
r 3 

2 ° 

5ab = r 1 AB , 
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and thus the non- vanishing connection coefficients are: 
1 



rr - |(S V Y r dr9rr 



(r - r c )(r + \r c \) 



1 1 

2 r (rc - r H ){r-H + \rc\) r-r H 

_ 11 rc (r - r H )(r + |rc|) 

2 r (r c - r w )(r c + |fc|) r - r c 

|fc| (r-r«)(r-r c ) 



1 1 

+ 77- 



2r (|rc| +r«)(|rc| + ^c) 
1 



r + r w 



1 1 1 

+ + 



1 

r r — r-u ' r — rc ' r + \rc\ 



and T[ t as well as 



A o 
\b = 3 (r - r w )(r - r c )(r + |r^|) Tab 

V 

-04 ■ 



Let us first consider 



We have 



rA r u A 



or 



(M')^rr _ fg- 1 y r Y r 



(M')tt 



l^^^ T ^t 

rt 



and 



Now, 



w/ B = ^( 5 -y B : 



= l/ % 1,3,2 ,&tp. 

rr 2^0r' 2 V (r - r H f(r - r c ) 2 (r + |f^|) 21 9i J 

1 3 r 
+ 77-r 



2 A (r - r«)(r - r c )(r + |r c | 



|W| 



- r n )(r - r c ){r + \r c \)\p^\ 2 



and 



A (r — r w )(r - r c )(r + j^ffj ' 7% ' 
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also note that 



1 ~ 1 1 A . . . w . „ / dib \ 2 

;T„ = ---(r - r n ){r - r c )(r + 



2r 3 



1 3 



r (dip ^2 1 1 / 1 2 

2A(r~r H )(r-r c )(r + |r^|)^ j 



1 T /dtps 2 1 r / 9 1p n 2 1 I — , 1 2 

2 L ^) + 2^ + 2> W I ' 



We then find that 

if o + Ki + K 2 



K M ' = I 



Kq + Ki- K 2 



dip ^2 1 



where 

JTn = 41 



(r - r c )(r + \r c \) 



2 r (r c - r H )(r n + \rc\) r - r u 

11 ~ 7^)0 + |re|) 

2 r (rc - r n )(r c + \rc\) r - r c 

_ 11 l^cl (r - r n )(r - r c ) 

2 r (|fc| + r H ){\rc\ + r c ) r + \rc\ 



1 /Ar 2 2m 



2V 3 r 



d 



i 



3 ' r 
1 3 



> Ar" 2m 



2 A (r - ru){r - rc)(r + \rc\) ' 



and 



and finally 



2r l 



r — ru r — rc r + \rc\ 
1 

#2 = 2-. 



Next we observe that 



2 r (r - r H )(r - r c )(r + \r c \) 
6 m 



2r 3 - ^ + (r - r H ){r - r c )(r + |fe|) 

- r(r 2 + \rc\r - r c r - r c \rc\) 

- r{r 2 + \rc\r - r«r - ru\rc\) 



r(r 2 - r c r - r H r + r H r c ) 



by (|3~TT|) , and thus 



/m; + Av i - - h~> = 2K t K, = - f-3 + — + — + —r=; 

r\ r — r-n r — rc r+\rc\' 
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1 1 

r (r — ru){r - rc)(r + \rc\ 
3 2 1 



ItuVcV + 2r c \r c \r - 2r H r c r - Zr c \rc\ru 
3 2 1 



A r (r — ru){r - rc){r + \rc\) 
again using the relations (|3.11[) ; and finally 



r — 3m 



> 



A r (r - ru){r + \rc\) ' 



K + Ki + K 2 > 2K 2 = 4- . 

r 



We have shown in particular 
Since 



K M ' > . 



r 2 yrr 1 r 



(M)^ t = I(M')^ 

r 

^ " AB = ~ y 'VAB , 

r 



we conclude 



K 



M 



1 



(.9 



-l\rr 



T ■ 



-K 



M' 



= ^\(r~ r n ){r - r c )(r + \rE\)T rr + -K M ' > ±±T r „ 



□ 

It is now an immediate consequence of the following (slightly more general) 
Gronwall inequality that (|4.5[) implies the weighted energy inequality (|4.6|) . 



Lemma 4.2 (Gronwall inequality for decreasing functions). Let a < 0, and f,g <E 

C 1 ([ri,r2]) with g 0, satisfying the inequality 



f < a y -f 
9 



(4.16) 

on the interval [r% , r 2 ] . Then 

( 4 - 17 ) /(r 3 )<{4^S/(ri). 



| 5 (r 2 )|H' 



Proof. By (|4~T5)) . 

(4.18) 



d 

dr 



/(r) exp 



^dr' 



< 0. 



which yields (|4.17p upon integration on the interval [ri,ra]. 

Finally, to verify that (I4.6P implies the statement of Prop. l2.1l note here 



□ 



1 1 



(4.19) J M -n=--T rr = 

r <p 



1 1 ,/A 



2r 2 



2to — r 



dtps 2 
dr ' 

1 ,dip,2 11 

1 + 27' 



2^4|i + 2 m-r V 9^ 
which completes our proof of the energy estimate in the expanding region. 
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Remark 4.3. The discussion of this section can be easily generalized to solutions of 
the Klein- Gordon equation 

(4.20) n g i/j = m A ^ (m A > 0) , 

by consistently replacing the standard energy current (|4.9|) by the modified current 

(4.21) J M - mA = J M - M b , 

where M b is simply the 1-form corresponding to the vectorfield M, i.e. ■ X = 
g[M,X). A detailed account of this generalization can found in |Schl2] for the 
example of the Klein-Gordon equation on de Sitter cosmologies; c.f. Section [51 

4.2. Local redshift effect. In this section we construct a vectorfield that captures 
the local redshift effect of the cosmological horizon. The basic insight for this 
construction was gained in the context of black hole event horizons in [DR08 . 

The idea is to construct a strictly timelike vectorfield TV in the vicinity of C + 
which is Lie transported by T, [T, N] — 0. As a consequence of the positivity of the 
surface gravity this can be done in such a way that the divergence of the current 
J N is positive. Thus, with E and the notation as in Prop. 12. 2[ the energy identity 
for J N on 

j+(£'uc+)nJ-(£ ro ) 

implies for ro > rc small enough (since the positive definiteness of the divergence 
only holds locally) 

(4.22) f J N -ndfig r < [ *J N + [ J N -ndng, 
Js ro n3+(s) '° Jc+ Js 1 



which leads to Prop. 12.21 given the left hand side provides an estimate for the 
boundary terms arising from the current J M of the previous section. 

We shall now define N in the double null coordinates of Section 13.21 On C + 
we may introduce the null frame (T, Y, : A = 1, 2), made up of the vectornclds 
T p.30|) . Y\ c+ as discussed in Section I3T41 and complemented by an orthonormal 
frame field Ea : A = 1, 2 on the sections of the horizon which is Lie transported by 
T. We have seen that, also using p.29|) . 

(4.23a) T = Kc u-^-, 

c+ ou 

d Id .41 



— L c K c guv — —2 , 

c+ 



(4 ' 23b) y c+ - JF^^- ic ^< ic -- c Wc + 

are conjugate null vectors on the cosmological horizon 
(4.24a) g(T,Y) 

(4.24b) g(T,E A )\ c+ =0 g(Y,E A ) c+ =0, 

and satisfy the commutation relation (|3.38l) . 

Let us now extend the vectorfield Y away from the horizon by 

(4.25) VyY = —(t(Y + T) , 

where u > 0. It is well known from [DR08] that this extension gives rise to a 
positive current by virtue of Lemma 13.11 Indeed, we find using p.37[) and (|4.25[) 
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that 

(4.26) K Y 



+ \°W\ 2 + l{T-^)(Y-i>) 



and thus 



(4.27) 



K 



> \kc{Y ■ ^) 2 + \\a (-) 2 -l (T ■ rpf + Un\ 2 > • 
C+ 4 v ' L2 y r c n c i 2 1 1 



if er > (2/rc) 2 2/«c. By construction Y is Lie transported by T, 



(4.28) 



[T, Y] = . 



Proposition 4.4. Let E, E', and &e as m Prop. I^.^l and let tp be a solution 
to (|2.7j) . Set N = T + Y with Y as constructed above, then there exists ro > rc 
(depending on A, m and EJ such that 

(4.29) K N [ip]>0 on J + (E'uC + )nr(E ro ), 
and a constant C that only depends on ro such that 

(4.30) J M [ip] -n< C(r ) J N [tp] ■ n on E ro nJ + (E). 

The positivity of the divergence (|4.29[) follows by continuity from (|4.27p with 
ro — rc chosen small enough, and a as indicated above. Also by continuity N 
remains timelike in a neighborhood of C + for ro — rc > small enough, which 
implies (|4.30|) : note that the uniformity of the constant C in r follows from the 
invariance of N under the flow generated by T, namely (|4.28[) . 

Alternatively, we can solve the O.D.E. (I4.25[) with initial conditions (|4.23[) in the 
double null coordinates of Section 13.21 If we set a 



2( — ) — , we obtain 

V rc I kc ' 



(4.31) N = k c 



2,2,3 
— ( — J uv 
kc rc 



d 



ic 



1 

uv 



K£_ 



du 

lc rc k c 



„ 1 A , o 

2 H (rl 

. 4 3 1 c 



ru\rc\) 



which allows us to verify all the statements of Prop. 14.41 explicitly; see [Schl2] for 
details. 

This completes the proof of Proposition 12. 21 



4.3. Limiting quantity. While the bound (|2.12[) follows from the global and local 
redshift estimates as discussed in Section [231 in particular (I2.10[) , we remark that 
the existence of the limit follows from a density argument. Any data ip\s € H (E) 
can be approximated by smooth functions tp . € C°°(E). Using the Sobolev inequal- 
ity on the sphere at infinity (Lemma lA.2[) . and commutations with the vectorfields 
f2(i) (|A.4[) . we can show that all solutions ipj to (|2.7p with initial data ip on E 
have a limit on E + , or on uv = 1 in the topology of M 2 of the Penrose diagram in 
figure[2l by integrating from E to E + as discussed in Section l4~5l The trace of tp on 
E + is then approximated by the limit of ipj in the norm defined by (|2.12p . namely 



the homogeneous H 1 norm on the cylinder 
Theorem O 



This completes the proof of 
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Figure 5. Localisation of the global and local redshift effect. 



4.4. Localisation. In this section we show that quantitative energy decay along 
the cosmological horizon is propagated through the redshift region and inherited 
on the future boundary. 

The observation that quantitative energy decay is preserved in a "redshift region" 
was first made by Dafermos in the context of the nonlinear Einstein-Maxwell scalar 
field model in spherical symmetry |Daf05j . It was applied by Luk to the study of 
the wave equation in the language of vectorfield multipliers [LuklOj . who showed 
that any polynomial decay rate in the exterior of the Schwarzschild black hole is 
inherited on the event horizon (and in a neighborhood of the horizon in the interior 
of the black hole); it also appears in Dyatlov's argument |Dyall| on exponential 
decay for solutions to the wave equation on Kerr de Sitter. 

The global redshift effect that is present in expanding spacetimes extends this 
mechanism beyond the vicinity of the cosmological horizon. We give here a proof 
which avoids a recourse to a bootstrap argument, previously used in [Luk 10 . 

Let us first look at the local redshift effect near the cosmological horizon with 
arbitrary incoming energy flux. Recall the notation introduced before Cor. l2.9l and 
the vectorfield N of Prop. 14.41 

Proposition 4.5. Given ip in the class of finite energy solutions to (|2.7[) as dis- 
cussed in Thm. \2.7\ let us define an incoming energy flux density g by 

(4.32) f *J JV M= r 9 (r)dr, 

where T2 = \Ct 2 > f or an V r 2 > T i- (Note that the left hand side is finite for 
any t-i > t\ by Prop. \2.S\ ) Let tq > rc be chosen according to Prop. \4-4\ and define 



(4.33) f(r) = / *J N &} 



c; 

where C' T = C r l~l J~ (E ro ) denotes the segment of the outgoing null hypersurface C T 
in the past of S ro (c.f. figure^). Then 

(4.34) /(r 2 ) < /(n) e - 6 ^- ri ) + f * g(r)e- b ^-^dr , 

Jti 

for all T2 > Ti, where b > is a constant that only depends on A. and m. 
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Remark 4.6. If we have no incoming energy flux from the static region, g = 0, then 
the energy / decays exponentially and proportionally to the strength of the redshift 
effect as reflected in the constant b > 0. In general, g > 0, the decay of / cannot 
be stronger than that of the incoming energy flux. Note here however, that the 
second term in the upper bound is a convolution around r = r 2 , and thus localises 
the contribution from the incoming energy. 

T N 



Proof. Consider the energy identity for J on the domain bounded in the past by 
C+ U C Tl and in the future by C T2 and E ro , r 2 > t\ , namely 



*J N + / J N -nd^ + K N dn g < *J N + *J N , 

where E^ ,Ta) = S ro n J + (C Tl ) n J~(C n ), and lli^X = J+(C+ UC Tl )n J _ (C^ U 
Ero 1 ' 1 " 2 )■ Observe that in the double null coordinates of Section I3T2I 

d/i s = — r dn A dt; A d/ia , 
where f2 2 on C + is a constant that merely depends on A and m, 

91 4 5 



Moreover, in view of (|4.23|) and using (|4.27j) . 

if L>mm -K C t c ,2 - - ' 7T c+ ■ 

Therefore there exists a constant b that only depends on A, and m, such that 
/ K N dfi g > b(A, m) j 2 dr / * J N , 

where we have used that by definition 

dr, 11 

du lc+ kc it 
Thus the energy identity takes the form 

f(r 2 ) + b f 2 f(r) dr < C g(r) dr + /(n) , (r 2 > n ) . 



This is a Gronwall-type inequality which easily implies (I4.34[) . Indeed, if we set 

F(r 2 , n )= / T2 {6/(r)- 3 (r)}dr, 

then the inequality simply reads 

/(r 2 )+F(r 2 ,r 1 ) </(n). 

Furthermore, 
d r 



9r 

so by integration, 



F(r 2 ,r)e- b ( T2 - T ) > -&[-/(r 2 )] e -^- T ) - g(r)e- b ^- r) 



F(T 2 ,T 1 )e- b ^ 2 - Tl '> > /(r 2 )(l - e^ 2 " 7 " 1 )) - / . 9 (T)e- b ( T2 - T W 
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which yields (|4.34p after rearranging and using the above inequality once again. □ 
Let us now complete the proof of Corollary 12.91 



Proof of Cor. \2.9l In a first step we apply our Prop. l4~5l to the local redshift region. 
Under the assumption (|2.13|) we are in the situation that the incoming energy flux 
density is bounded from above by 

9(r) = ^ (fc€N), 
and similary under the assumption of exponential decay by 

g (r) = fie-** (/? > 0) . 

Now we have, on one hand, 

1 Jt i Jti ' 

< — + e~ aT2 — = — + e~ aT2 — 

where we have chosen t[ £ (t~i,t 2 ) such that r 2 — r{ = \ar 2 for some a £ (0,6), 
and on the other hand, 

f irV"^ 2 >P< b 

pe-^e-^-^dT < I /3(t 2 - Tx)e- bT2 , = b . 

We conclude by Prop. 1431 if (|2.13p holds then also 

/ *J N m < % + e- b ^-^ f *J N [>p} (r 2 > n) , 

JC' T2 T 2 JC' ri 

and similarly if (|2.15l) holds, then there exists a constant C (we assume for simplicity 
(3 6) such that also 

*J N W\ < Ce~ min Wrt T2 +e - b{T2 - Tl) f *J N [iP] (r 2 > ti). 

In a second step we apply our Prop. 14. II to the global redshift region. Here we 
localise the argument of Section 14.11 to the domain of dependence of C+ U C' T (as 
indicated in figure [5]). Let ro be chosen as above, and denote for r > ro by 

K = (s r n J+(c T )) \ j+(c t \c;) 

the segment of E r in the domain of dependence of U C^. for any given t > tq. 
Then, as discussed in Section 14.11 c.f. in particular 



r 2 / J M -ndn-g < n / J -ndfig (r 2 > r x > r ) ; 

note here that the vectorfield M is timelike and hence gives rise to a positive 
boundary term on null segment of the future boudary of the domain of dependence 
of C+ UCj, which preserves the inequality (|4.3j) with XV replaced by E£.. Next we 
choose n = ro and consider the energy identity for J N on the domain bounded in 
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the past by C+ U C' T and in the future by E ro , as in the proof of Prop. 14.51 By 
Prop. l4~4l we then have 

r f J M -nd^ r <C [ *J N , 
J K Jc+uc T 

which decays according to our assumption and what is shown above. We may 
now take the limit r^ — > oo as discussed in Section 12.31 to obtain the statement 
of the Corollary with r replaced by r', where £+, = (£+ n J+(C T )) \ J+(C T \ 
C' r ). It remains to show that r' and r are related linearly, which is an immediate 
consequence of geometric relations described in Section 13.21 For given values of u 
and ro the incoming null hypersurface originating from the sphere at (u, tq) is 

_ho , _ r o - r c 

(rO - r H ) r C r H + \r c \ ( rQ _|_ | rc |) r c lr c \+r n 

which terminates on S + at the sphere (u',v), where 

, 1 1 
u = — = — u . 

v h 

Therefore, since T ■ r = 1, also 

1 f u ' 1 1 1 
t' = tH / — du = t H log — . □ 

n c J u u k c h 

4.5. Pointwise estimates. The usefulness of the foliation (E r ) is revealed in the 
Sobolev inequality for E r which allows us to pass immediately from energy estimates 
to pointwise estimates using the symmetries of the spacetimes (A4(A,-m)> <?)• 

Proposition 4.7 (Sobolev inequality on E r ). Let ip £ H 3 (E r ), r > rc, then 



Ar 2 2m 9 , 9 

1 1 r sup i/r < 

3 r E 



< C 



3 3 

/{^ 2 + E(%)V0 2 +E(%)%)^ 



3 

, 2 \ - /_ _ , s 2 



+ (tv>) +^(%)T^) + ^ (n w n (i) TvO |d^.. 

on any submanifold E C E r o/ </ie /orm / x SO(3), where I CM. is an open interval 
with \I\>1. 

The angular derivatives on the sphere are here estimated using the generators 
of the spherical isometries Qr{\ : i — 1,2,3. This fact is known as the coercivity 
inequality on the sphere which is discussed in Appendix \K\ where also the precise 
definition of the vectorfields fi(j) can be found. 



Proof. Recall the induced metric (|2.4I) on E r . Let t , to £ I, and £ £ § 2 , then 

dtp 
~dt 



2(5 
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Therefore by Corollary IA.3| there exists a constant C > such that 
|rVW-rV(to)|< 

< 



"' I "' S r i=l i,i=l 

PC 3 2 3 2 



+ c 

Now choose (to, € / x S 2 such that 

and we obtain the stated inequality in view of the expression (|2.5|) for the volume 
form of g r . □ 

The application of Prop. 14.71 to the functions 4>%T ■ ip and <f>^ ^tt^ip ■ i = 
1,2,3 yields quantities on the right hand side which are monotone by our results of 
Section |4~T1 provided ip is a solution to the wave equation (12 . T|) . Indeed, we have 



(4.35) sup<^ (r</>) 2 {T^f + r 2 \fib\ 2 \ < 



< C 



f {^{T^f + 1^| 2 + ]>> 2 (™ w ^) 2 + Y>o w ^| 



3 

, 2 



3 3 n 

+ £ </» 2 (T 2 fi (i) (j) ^) 2 + 51 |^%)^(i)^r Ud/i Jr , 

which is bounded by Prop. l2~Tl because T and Q(j) are Killing vectorfields, and thus 
£l(i}ip, Q,n\Clrj\tp, and Ti/j, TQ,n\i/j, TQ,n\Q(j-\iJ) are solutions to the wave equation. 
In view of 



3 > 
lim (r<f>) 2 = — , r 2 |)(7^| 2 = IX 7 V| 



the pointwise estimate (|2.21[) is now obtained analogously to Thm. 12.71 as discussed 
in Section l2~3l 

We may also apply Prop. l4~7l to ^§7r</> 3 which yields by Prop. |2~T1 a monotone 
quantity. In fact we obtain, 

(4.36) (r(/ ,)-2 sup (^) 2 <^ cMi 

s r or r* 

which implies by integration from E along i-constant lines the final bound (|2.20p . 
This completes the proof of Corollary 12. Ill 
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Figure 6. Penrose diagram of de Sitter depicted as a spherically 
symmetric spacetime with respect to a chosen timelike geodesic T. 

5. Epilogue: The wave equation on de Sitter spacetimes 

The de Sitter cosmology is homogeneous and can be cast as a hyperboloid in R 5 
with spatially compact sections S 3 . The global study of solutions to the wave 
equation reduces in this setting in principle to a local problem due to the presence 
of additional symmetries. In fact, it suffices to localise to the past of a point on the 
future boundary to obtain a global result; see also the work of Vasy |Vas07j and 
[Bas09|lBasl2] . 

In view of applications to Kerr de Sitter cosmologies however, we may entertain 
an approach to the global study of linear waves on de Sitter which does not make 
use of the homogeneity of the spacetime. This consists in viewing the spacetime 
as spherically symmetric with respect to a fixed timelike geodesic, and corresponds 
precisely to the case m = (see Section l37T| or |GH77] ). It does have the advantage 
that a global causal geometry emerges that shares some essential features with the 
Schwarzschild de Sitter geometry, as depicted in figure |5] in particular, there exists 
an expanding region in the sense described above, where the metric takes the form 



3' 

which is bounded in the past by cosmological horizons, and to the future by a 
spacelike hypersurface E + at infinity. 

This approach is carried out in all detail in my thesis jSch!2j . The analysis 
yields results analogous to those of SectionEJ and can form the basis of perturbative 
treatments. In particular we establish Theorem 12.71 in the case m = 0. 

Moreover, in Sc hl2j we provide an alternative proof of Prop. [2TT1 in double null 
coordinates, using the global redshift multiplier 



which highlights the global redshift effect as a suitable extension of the local redshift 
effect, (c.f. Section l4~2j). 

Finally we note that our approach lends itself to the study of global solutions to 
the Klein-Gordon equation O g ip = m Aip, tua > 0, by introducing consistently the 



(5.1) 



g = -<p 2 dr 2 + g r 



1 



dr 2 +g :i 



(5.2) 
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modified energy current 




we prove in addition to Thm. 12.71 



that the limit on S + vanishes identically as a L 2 (R x § 2 ) function |Schl2j . 



This paper was completed as part of my thesis. I would like to thank foremost 
my advisor Mihalis Dafermos for many insightful discussions and for suggesting 
this problem to me. I am very grateful for his support during my Ph.D. I have 
also benefited greatly from conversations with Willie Wong and Jared Speck during 
their time at Cambridge, and enjoy the frequent exchange with Stefanos Aretakis. 
The author is supported by the UK Engineering and Physical Sciences Research 
Council, the Cambridge European Trust and the European Research Council. 



In this appendix we shall prove the coercivity formula for the standard sphere, 
and recall the classic Sobolev inequality on the sphere. 
Let us denote by 
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Appendix A. Coercivity inequality on the sphere 



(A.1) 




the sphere of radius r in Euclidean space, a submanifold of 



(A.2) 




We denote the metric of the round sphere S r by 
(A.3) 7r = e | T ^ =r 2 7, 

where 7 is the standard metric on the unit sphere § 2 . 



Let 



(A.4) n (i) = g. ifc a?J i = l,2,3, 



where e is the volume form of e. We have 



3 



(A.5) 




and thus for all x £ R 3 \{0} and u : R 3 -> 



E 



diffcrcntiable: 



3 



(A.6) ^(% )U ) 2 (x) = M 2 [|Vu| 2 - (iL v u > 2 ] = |*| W . 
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Lemma A.l (Coercivity inequalities on the sphere). Let u be a smooth function 
on S r , then 

(A.8) r 2 |JHl<I>(^) 2 , 

i=l 

and 

3 

(A.9) AV 2 u\l< J2(%)^U) U ) 2 - 

»>i=i 

We have already shown the first inequality. For the second inequality we can use 
that more generally for any SV-l-form (a 1-form on R 3 such that 8 ■ X = 6 ■ HX) it 
holds 



(A.10) EK%>C = 



3 

2 21 

r 



where denotes the Lie derivative on S r . (To prove (|A.10[) one can proceed analo- 
geously to Lemma 11.2 in [Chr09 .) By substituting 

(A.ll) 6 = iu = Au\ TSr , 

we then obtain 

(A.12) r 2 \p 2 u\ 2 =r 2 \piu\ 2 < 



3 3 3 

< El^ o w Ht = El^ n (0 < - EIF(%)l 7r - 

i— 1 i— 1 z— 1 

because Lie derivatives commute with exterior derivatives. Inequality (|A.9[) then 
follows from (|A.12I) using (|A.8|) . 

We recall the classical Sobolev inequality on the sphere. 

Lemma A.2 (Sobolev embedding on § 2 ). Let u £ H 2 (§ 2 ), then u G L°°(§ 2 ) and 

(A.13) ||u|| L oc( S 2 ) < C||U|| H 2(S2) ■ 

Given a function u on S r we can apply Lemma lA. 21 to uo h, where 
(A.14) h:S 2 ~>S r , £H->r£. 

Since, in the coordinates 

(A.15) 7r = r 2 7= r 2 1 AB dy A dy B , 

h is the identity mapping, we get 

(A.16) \^7{uoh)\l^r 2 \^7u)oh\ 2 lT , & (u o h)\ 2 = r 4 \W 2 u) o h\ 2 T , 

and thus 

(A.17) r|t»|| s < C (J \ u \ 2d ^,y +C(J r 2 \^u\ 2 dfi Jr 



+ C(J r A \p 2 u\ 2 A^ lr 



30 



VOLKER SCHLUE 



Corollary A. 3. Let u G R 2 (S r ), then 
(A.18) r\u 
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